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1 Introduction 



Interactions of pions at low energies are dominated by the chiral symmetry of 
QCD. The best framework for this is Chiral Perturbation Theory as worked 
out systematically for the two-flavour case in JT|. The assumptions involved 
here are that the Goldstone Bosons, identified with the pions, resulting from 
the spontaneous breakdown of chiral symmetry are the only relevant degrees of 
freedom. We can then expand in the small quark mass and small energies and 
momenta. This expansion is good up to a scale of about the chiral symmetry 
breaking scale of order of the p mass. This expansion works very well in the 
purely mesonic sector, some reviews can be found in [0. 

In the present article we will deal mainly with the SU(2) ® SU(2) chiral 
symmetry, where the isospin limit m u = is taken. For this case there exist 
already several full two-loop calculations 0, |], |^| and also some calculating only 
the dispersive parts|| 0. The expansion parameters x ^ F i and |if are quite 
small but still it is sometimes useful to go beyond 0(p A ) because of the following 
reasons: 

1. The results at 0(p A ) do not fit the experimental data like in the 77 — > 7r°7r° 
case@. 

2. One deal with pions in the final state. In this case a isospin zero S-final state 
produce a strong rescattering, and one can obtain huge corrections [|. || . 

3. The quantity under consideration itself only starts at 0{p A ) so in order to 
have an estimate of the accuracy it is necessary to go to 0(p 6 ). 

The calculation here is case 1, see section pT2] , and 3. 

There are two problems involved in going beyond 0(p A ): 

1. The number of constants increases rapidly with the order. At 0(p 2 ) the 
lagrangian contains two constants, at 0(p A ) 7 and at 0(p 6 ) there is only a 
classification for the three flavour case||. 

2. The algebra involved in such a calculation is rather tedious. 

Problem 1 we address by using resonance saturation estimates of the relevant 
constants |J. This is known to work well at 0(p A ) and was also used in several 
other two-loop calculations. 

This paper is organized as follows. In Sect. ^| we define the form factors and 
discuss the present experimental status. The next section reviews the presently 
known theoretical results, which we have checked independently, about A and 
V. Section |4] presents our main result, the two- loop calculation of A. Technical 
aspects have been placed in two appendices. We then present numerical results 
and our conclusions. 
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2 Definitions and present experimental results 

A review of the theoretical and experimental situation up to 1982 can be found 
in 0. 

2.1 Matrix element 

We will consider the ir~ decay 

7T~(p) -»• l~(pi)v{Pv)l{q) [^2 7 ] (!) 

where / stands for either e or /x and we will deal with a real photon, i.e (q 2 = 0). 
The process tt + — > Z + z/7 can be obtained by charge conjugation. The decay with 
a muon in the final state is completely dominated by Bremsstrahlung, while this 
contribution is helicity suppressed in the electronic mode. This decay therefore 
presents a good place to look for the pionic structure. 
The decay width can be written as: 

dT = 2M (2ir)5 ^ \T 2 \d LIPS (p; Pl ,p u ,q). (2) 

n V / spins 

with 

dLIPsfaPl, ■■■,Pn) = $ 4 (P - ^2 Pi) II 7?) ( 3 ) 

i=l i=l Z Pi 

the phase space volume, and we have used the covariant normalization of one- 
particle states: 

< p'\p>= (2n) 3 2p°5 3 (p' - p), (4) 
where the matrix element T is calculated inside the V — A theory [11 1 : 



T = -iG F eV u * d e;{F n L» - H^l v } (5) 

with: 

v = m ^)(i+ T5 )(^-M±iz! )ufe) 

2p ■ q 2pi ■ q 

F = u(p u )Y(l-l 5 HPi) 

= iV(W 2 )e^q a pp- A{W 2 ){q-Wg^ u -W^q u ) 
W » = (p-qY = (pi+p„Y. (6) 

Here e M stand for the photon polarization vector with g M e M = 0, whereas A and V 
are the most general two Lorentz invariant amplitudes occurring in the decom- 
position of the tensor : 

//- = J dxe n-*+iW-y K o\TV^ m (x)IU 5 (y)W + (p) >,I = V, A. (7) 



at q 2 = 0. The most general decomposition can be found in [12 
As usual we divide the amplitude in two pieces: 
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1. The inner bremsstrahlung (IB), the photon is emitted by an external par- 
ticle. The pion or the lepton in our case. This corresponds to the term 
containing F n L^ in (|^). This is the radiation of a pointlike electron and 
pion. 

2. The structure dependent part (SD), it is in these terms that the bound 
state structure of the pion plays a role. This contains the two structure 
functions V(W 2 ) and A{W 2 ). 

Usually the internal Bremsstrahlung plays the mayor role, masking the struc- 
ture dependent effects. But in n — > ewy it is helicity suppressed, allowing de- 
tection of the structure dependent terms. A more extended discussion of the 



kinematics can be found in the review of particle properties |13fl or Ref. [p^| . 

Time-reversal invariance implies that A and V are real functions for W 2 below 
the two-pion threshold, which is the region of interest here. They are analytic 
functions of W 2 with cuts on the positive real axis. 

One of the reasons to perform the present calculation is that the W 2 depen- 
dence of the two form factors only starts at 0(p 6 ), see below. 



2.2 Present experimental results 

There are two experiments that have determined ^(W 2 ). They both assumed a 
constant form factor. One of the experiments |14| used the decay tt + — > e + ve + e~ 
thus allowing to determine the sign of V as well. The average quoted in |T3| is 

F v = (0.017 ±0.008). (8) 

The axial form factor is only measured via the ratio 7 = A/V. Assuming the 
CVC value F v = 0.0259 ± 0.0005 this yields@ 

F A = (0.0116 ±0.0016). (9) 

These related to the ones we use by 

-V2m n+ (V,A) = (F V} F A ) . (10) 

There are in the literature some contradictory experimental results involving 
the pion decay in a photon and a semileptonic pair. The last one was reported 
by Bolotov [K] .There they explore a wide kinematical region, observing 80 7r 
decays in flight, and obtain a result that for the decay width deviates from the V-A 
standard model theoretical calculations by more than three standard deviations. 
This results disagree strongly with previous ones [0 , that are in agreement with 
the theoretical calculations (keeping uncertainties). To obtain a good fit to the 
data in ]T3] they introduce a tensor radiation term of the form: 

Ttensorial = f ^ * & V ' F T u{p e )o >(l + l 5 )v(Pv) (H) 
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with a coupling constant Ft = —(5.6 ± 1.7) x 10~ 3 . This terms interferes de- 
structively with the usual V — A form. 

The tensor force would be a signal of new physics. A lot of effort has been 
devoted to explaining this experiment, see e.g. [|17]], where supersymmetric ex- 
tensions of the standard model have been used. This mechanism is not able to 
explain the huge disagreement between the data and the theoretical calculations. 
One can try to learn something from similar decays more accessible experimen- 
tally, see e.g. [18]. Looking at K + — > 7i ue + as was done in ]19| could not rule 



out the possibility of a small tensorial term in the amplitude. 

This motivated us to check whether there were any anomalously large mo- 
mentum dependent effects in the form factors. These are not expected [|TTJ but 
a definite calculation of this effect within CHPT was not performed up to now. 



3 Previous CHPT results 
3.1 CHPT Lagrangian to 0(p 4 ) 

We use in this work the sigma model parametrization as used in 0J. This uses 
the equivalence of the groups SU{2) ® SU (2) and 0(4) and simplifies the vertices 
compared to the exponential parametrization. 

The most general effective lagrangian consistent with parity, Lorentz invari- 
ance and chiral symmetry is given at lowest order by: 

F 2 

C 2 = —V^V^U + 2BF 2 (s°U° + p l U l ). (12) 

Where B and F are constants not fixed by symmetry. We will work here in the 
Standard CHPT assuming a large value of B. U(x) denotes a four-component 
real field in the vector representation of 0(4). s(x) and p(x) are the scalar and 
pseudoscalar external fields respectively. We introduce the external sources in 
the usual way pi] . The covariant derivative is defined as: 



Vy7° = d^U° + a^{x)U l 



\/^U l = d^U 1 + e ikl v k {x)U l - a? (x)U°. 



Here v h (x) and a l (x) are the external vector and axial- vector sources. 
At 0(p A ) the Lagrangian is given by 



(13) 



£ 4 = ^(V^V^f + ^V^V^XV^V^) 
+h{x ] Uf + Z 4 (vVV„E7) + h^F^F^U) 
+h{V^F, u V u U) + hix^Uf + h lX *X + h 2 F, v F^ 
+h 3 X ] X- (14) 
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Where the tensor F^ u is defined by: 

(V M V„ - V V V,)U = F, V U (15) 

which contains the external fields v^, a M together with their derivatives. Also we 
have introduced the vectors: 

X A = 2B(s°,p*) 

X A = 2B(p°,-s l ) (16) 

We have used dimensional regularization in the MS scheme. In this regularization 
the low-energy constants are defined as: 

k = ll+jiX (i = l,...,7) 

A = l?^-5" og(to) + r ' (1) + 1)} (17) 

where the l\ are the coupling constants renormalized at the scale /i, uj = (d— 4)/2 
and the 7« factors are found via the Heat-Kernel expansion and are given in [l|]. 
The derivative of the V function is the Euler constant, T'(l) = —7 and later we 
will also use m 2 for the pion mass squared, m 2 . For the quantities considered 
here it makes no difference to 0(p 6 ) if we use the full pion mass or only the first 
term in its quark mass expansion. 

3.2 A form factor at 0(p 4 ) 

This calculation was performed in Ref. |l|. The loop contributions, though 
allowed by power counting, cancelled and the counterterms were the only nonzero 
contribution: 

A=j(2l & -Ah) . (18) 

This is in fact the experimental data used to determine this combination of pa- 
rameters. It is therefore also useful to have an estimate of higher order effects to 
know the accuracy of this determination. 

3.3 V form factor at 0(p 6 ) 

This part start to play at 0(p A ) and it involves at least one vertex given by the 
Wess-Zumino (WZ) lagrangian [KJ]. So due to the initial order at most one- 
loop diagrams can contribute to 0(p 6 ) constructed with one vertex coming from 
the WZ Lagrangian and the others from the lowest order chiral lagrangian (see 
below). We performed the calculation in the two-flavour case and our results 
agree exactly with those obtained by Ametller et al.|HJ] when restricted to two 
flavours. 
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The VfW 2 ) form factor obtained at 0(p 6 ) is: 



V(W 2 




m log( 



,771 



Am 2 
~W 2 ~ 



log( 



1 - 



4m 2 
H/2 



1 - 



4m 2 



- 1 



3 9 



+ 



+ 



(19) 

to this process 



The Kaons are frozen out and the full kaon contribution to 0(p 6 
can be absorbed in Z 4 , which in turn is absorbed in F n . 

The divergent parts also agree with those obtained in [22] for Np = 2. The 
finite counterterms were estimated by resonance exchange. We use the "hidden 
symmetry" scheme [23] for the vectors and calculate the possible contributions 
assuming full vector meson dominance. We have disregarded the contributions 
of heavier resonances than the p. 

Numerical results concerning the ^(W^ 2 ) form factor can be found in Sect. [5]. 



4 A form factor at G(p ) 
4.1 Overview 

All quantities have been computed using a Feynman Diagram technique. We 
have used two techniques, a brute force one where everything was calculated and 
the formfactor A extracted afterwards. This has provided an independent check 
of the values for F n and M 2 given in |J. The other technique was a more aimed 
extraction of the formfactor A. This is the one we will describe below. The A 
formfactor is the only one that has a contribution proportional to g^p ■ q. 

The presence of the g^ v requires that the axial- vector insertion and the vector- 
insertion are in the same one-particle irreducible subdiagram. This immediately 
removes a large part of the diagrams. The presence of the p ■ q kinematical factor 
then guarantees it is not part of the internal Bremsstrahlung contribution. 

The direct calculation of the A formfactor has the additional advantage that 
some of the most difficult aspects of renormalization at two-loops do not appear 
since the 0(p 4 ) contribution only was a finite counterterm contribution. We 
therefore only need the pion wave function renormalization, F v and m 2 to 0{p i ) 
accuracy. This is similar to the situation in the calculation of 77 — > 7r 7r° of 0. 

Another simplifying fact is that in the sigma model parametrization used here 
there is only a vector-two-meson vertex in the 0(p 2 ) Lagrangian. 

A quantity needed in the remainder is 

(^ e » - T^i { 1 1 + " < - ^ + -< ~ + - log (?) } • < 20 > 
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Let us now discuss the contributions from the remaining diagrams at 0(p 6 ): 

1. Tree level diagrams with one 0(p 6 ) vertex. These are needed to perform 
renormalization and the finite parts we estimate using resonance exchange. 
This is described in subsection 



2. Tree level diagrams with two 0(p 4 ) vertices. These can be combined with 
the one particle reducible diagrams containing two separate one-loop sub- 
diagrams and those with a one-loop separated from one 0(p A ) vertex. To- 
gether they combine to applying wave function renormalization to the 0(p A ) 
result. They contribute 

^-yi'-FiR*) <2i) 

3. The two- loop diagrams with nonover lapping loops, shown in Fig. |]a and 
the one loop diagrams with a 0(p 4 ) vertex insertion on a propagator in the 
loop never produce a factor p ■ q and hence do not contribute to A. 

4. One-loop diagrams with a 0(p A ) vertex. These are in fact the main contri- 
bution numerically of the loop diagrams, see section [|. Analytically they 
contribute, 

I Tfl 2 

A=— (-20/ 5 + Wl 6 - 8h + Al 2 ) —— B (22) 
F 6 [a — 2) 

5. The pure two- loop diagram of Fig. [L|b can never produce g^ u p ■ q and does 
not contribute to A. 

6. This leaves now only the diagram in Fig. [TJc. Its evaluation is the most 
difficult part of this calculation. Its treatment is sketched in App. |A]. 

7. The 0(p 4 ) part contains 1/F, we change this to 

8. The contributions from the counterterms are split into a finite and an infi- 
nite part. The finite part we estimate using resonance saturation and the 
infinite part will cancel all the 1/u divergent terms. In addition we choose 
to renormalize such that all terms log(47r) — 7 + 1 disappear as well. 

Putting all of these contributions together we obtain: 
A(W 2 ) 

F,|^i<*-M 1 + i^ m,+ (J^*) 



1 



2 

m 



-20k + 10l G -8h + Al 2 ) — — B 

(d-2) 
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(a) 




(b) (c) 

Figure 1: One particle irreducible two- loop diagrams for A. The circle-cross 
means a vector current, the black box an axial current while the cross is used for 
a vertex with both the vector and axial current. 
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(23) 



where £>o is defined in ( [20]) and /„, a a are axial- vector-meson couplings defined in 
Eq. (|26|) . We have explicitly shown the single poles, which have to be subtracted 
via counterterms. The renormalization scheme used corresponds to the one in 
Ref. IJ. We subtract in Eqs. fl23[) and (p0| ) the terms in [• • -} sub and replace I4 by 
The A( £: |) function used here is in Appendix We have obtained those 
only in numerical form. The integrals can be computed rather efficiently A full 
analytical evaluation might be possible similar to tctt scattering but we have not 
been able to do all the integrals analytically. 

We apply the renormalization group equation to the chiral lagrangian to ob- 
tain one check on our results [P^j . There are no one- loop contributions propor- 
tional to divergent combinations of the /j. Therefore there are no 1/u 2 poles in 
the final result. These cancel as required. 



4.2 Resonance estimates of the 0(p 6 ) parameters 

We evaluate the resonance contribution following ||. Invoking P and C invari- 
ance, the relevant lagrangian can be written as 



4= E ^Kin{R) + £>Int(R) \, 



R=V,A 



with the kinetic terms 



1 /if 2 
C Kin (R = V,A) = --{V x R Xl N u R^ - -fR,uRn 



(24) 



(25) 



where (C) means the trace over C. 

We have described the vector and the axial-vector mesons in terms of the anti- 
symmetric tensor fields V^ u and A^ u , where we have restricted ourselves to the 
octet fields and Mr is the corresponding mass in the chiral limit. 
The interactions read 



C 



Int 



2V2 



f A (A^) + ia A (A„{u v , /H) 



(26) 



All coupling constants are real and we have used 



R 



(27) 
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where the covariant derivative acts on the octet multiplets as 

d^R = d^R + [r„, R] (28) 
The connection is defined by 

T M = ]^{u ] \d^ - i(vp + ajj\u + u[dn - i{v^ - a M )]u + } (29) 

The vector-axial field has been defined by: 

m m = iu ] DJJu ] = ul (30) 

and the external fields-strength tensor has been introduced via 

f^=uF^v)±v)F^u (31) 

and are associated with the external left and right field sources. 

Keeping in mind that the resonance propagator cannot decrease the chiral 
counting, the resonance exchange between R, pseudoscalar mesons and external 
fields has to be 0(p 6 ). As input constants we have chosen in the a\ resonance 
exchange the following ones [ 25 ]: 

a A ~ -6.66 10~ 3 , f A ~ 0.080 (32) 

where they have been calculated at leading O(Nc)- With these values p5f quotes 
a partial width for the process uj — > 7r + 7r _ 7r of 7.3 MeV to be compared with 
the experimental value 7.5 ±0.1 MeV and a good value for the decay a\ — > 717. 
In this case both terms have contribution to the vertex ir~ — ► 701 while only the 
one coming with f A is found in the a\ W~ . 

There's no b\ contribution due to CP, neither any scalar or tensor contributes to 
this order because of spin-isospin. 



5 NUMERICAL RESULTS. CONCLUSIONS 

We comment briefly our results about ^(W^ 2 ) form factor. This is found to be in 
perfect agreement with the one found in |H]] . We plot in fig. |2] this form factor 
versus the lepton pair invariant mass squared with the photon on mass-shell. We 
find that inside the allowed kinematical region the variation is around 4.5%. 

For the evaluation of the A form factor we first use the central values of the 
l r quantities in [p], p6[| . These are related with the l r ones via: 



7» n , i f™ 2 
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Figure 2: V{W 2 ) form factor at 0{p A ) (short-dashed line),(9(p 6 ) (dashed line) 
and 0(p 6 ) plus VMD contribution (full line) vs. the semi-leptonic pair momenta. 



fJi 


m p 


0.6 GeV 


0.9 GeV 


and F — > F w 
0{p Q ) 1-vertex of £ 4 
0(p 6 ) pure two- loops 


-5.95 10" 2 
-0.22 10~ 2 
+ 1.03 10" 2 
+0.53 10~ 2 


-5.95 10" 2 
-0.24 10~ 2 
0.88 10~ 2 
0.42 10" 2 


-5.95 10~ 2 
-0.21 10~ 2 
1.19 10- 2 
0.59 10~ 2 


Total 


-4.62 10~ 2 


-4.89 10~ 2 


-4.44 10~ 2 



Table 1: Some contributions to the t4(PU 2 ) form factor at W 2 = for Z[(m p ) of 
fl3"3p. Units are GeV^ 1 . Column 2 is at our standard fi. The other columns show 
the variation with \x. 
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At the m p scale these are given by: 



l[ = -0.00540 V 2 = 0.00567 l\ = 0.00560 l r 5 = -0.00553 l r 6 = -0.01381 . 

(34) 

The other inputs are the charged pion mass and F n = 0.0924 GeV. The size 
of the various contributions at W 2 = can be found in Table [I]. The term 
proportional to A is very small, about 5 • lO^GeV" 1 . The estimate of the 
resonance contributions vanishes at W 2 = 0. The pure two-loop diagram provides 
about 1/3 and the one- loop diagrams with one vertex of £4 about 2/3. All 
together the total correction is about 25%. The remaining subtraction scale 
dependence can be judged from comparing the results of the three scales used in 
Table @. 

The combination of l r 5 and 1$ that appears here is determined from A but 
the variation due to the uncertainty on 2l± — l 2 is rather large. A better result 
is obtained by noting that the combination 2/ x — l 2 is directly obtainable in 
the isospin one, spin one hit channel. It is therefore more suitable to use the 
experimental value of a\ directly. If we use a\ = 0.038 ± 0.002 and the value of 
l\ given above and the 0(p 4 ) expression for a\ of [l[] we obtain the range 

2Zi -l 2 = -0.0172 ±0.0072. (35) 

The contribution proportional to 21 1 — I2 from this range is 

A 2h _ h = (-1.84 ± 0.77) 10" 2 GeV' 1 . (36) 



A Roy equation determination of these constants using all available irn data|27| 
with h = -1.70 ± 0.15 and l 2 = 5.0 or 

2Zi -h = -0.0141 ±0.0003. (37) 

This leads to the more restrictive range 

A 2h _ h = (-1.51 ± 0.04) 10" 2 GeV' 1 . (38) 

The W 2 dependence is rather small and within the range relevant for this decay 
extremely linear. 

W 2 

A W 2 = (-0.0009 - 0.00110)— (39) 

rat 

TV 

The first term is the two-loop contribution and the second term is the resonance 
estimate. The total contribution from the resonance estimate to the formfactor 
is rather small, below 2%. 

Taking the particle data book values A = —0.0588 ± 0.0081 GeV -1 and the 
value in fl37"D we obtain 

2/5 -l 6 = 0.00315 ±0.00030, (40) 
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Figure 3: The A{W 2 ) form factor at 0{j> 4 ) (short-dashed line) , 0(p 6 ) (dashed 
line) and 0(p 6 ) plus VMD (full line) vs. the semi-leptonic pair momenta. 

to be compared with the value at 0(p 4 ), 0.00275, quoted above. The quoted 
error is only the experimental uncertainty. Additional uncertainties are the \x 
dependence (about 5%), the fact that the other l\ are estimated using F n = 
93.3 MeV (a few %) and contributions from other resonances. Adding those in 
quadrature would increase the error in (f40p to about 0.00035. We have plotted 
the A form factor in Fig. |3| to show the small W 2 dependence. For completeness 
we also give the ratio of the two form factors, 7 = y in fig. f|, together with the 
present experimental data. 

The 0(p 6 ) corrections to the A form factor are sizable and diminish it in 
absolute value. The W 2 dependence is very small as expected and is dominated 
by the resonance estimates. The deviation from the V — A picture observed in 
JOJ can therefore not be explained by the small W 2 dependence observed here. 

The rather large correction found here for A also affects the prediction of A for 
the decay K — > ^7 ||12|| . As a first guess at the size of the correction we have also 
calculated the correction with M K and F K used instead of m n and F n , this results 
in a 20% downward correction. For a definite prediction for the formfactors in 
the Kaon decay we have to perform a full SU(3) calculation. 

In conclusion, we have performed a two-loop calculation in CHPT for the A 
form factor in tt — > Zz/7 and estimated the relevant new constants by resonance ex- 
change. The corrections are dominated by the loop effects and the pure two-loop 
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Figure 4: The 7 = £ factor at C(p 4 ) (dashed line), and C(p 6 )+VMD (full line) 
the semi-leptonic pair momenta. We also shown some experimental results 
15| and m 



vs. 
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diagram is sizeable. The total correction is about 25%. We have also confirmed 
the known results for the V formfactor. 
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A Appendix A 

In this appendix we collect the main functions that enters in the evaluation of 
the A form factor at 0(p 6 ). 

We will take advantage of the fact that we want to extract the part propor- 
tional to g^ u p ■ q. With the notation defined in App. [B], the presence of g^ u 
requires the presence of r lAt and a factor of r lz , or r 2u from the other vertex. We 
then rewrite r\ ■ r 2 , r\ and r\ in terms of the possible denominators to put the 
most divergent parts in terms of integrals with two denominators only. Those 
never contribute a factor p ■ q. 

The remaining integrals are of five types: ((r lf ,r lu )), ((r lf ,r 2u )), ((n^i^ia)), 
(( r iti r iu r 2a)) an d {{ri^r 2v r 2a )) . The second and the fourth are trivially related 
to the first and third respectively by using identities on the r 2 subintegral. The 
procedure of evaluation is described in more detail in App. [FJ. 

The function Af^) is defined as follows: 



• rn 

A(^—^-) = ^[ [ [ dxi dx 2 x 2 x 1 (l - x 2 )\x 2 (-A + 3x 2 + 4x 

m z 6 J Am 2 a Jo Jo I 

+^-x 2 U + 6x 2 - 2Ax 2 2 + 2! (7 + 10x 2 - lOx^; 



+ K — f x 2 2 (-2A + 32x 2 ) + xix 2 (36 - 88x 2 + Uxj) 



(p • lY 

m 4 

+^(2-x 2 ))}f 3 [z] (41) 

As one can see this function make use of the F 3 [z] (defined in Appendix 
B), which is finite. We evaluate this functions numerically. A( £ 4) is a very 
slowly convergent function, so one must pay attention in its evaluation. We 
have first made a conformal transformation over the a variable, which makes the 
integration over the new variable simpler. The remaining two integrations (x±, 
x 2 ) have been evaluated in three different ways, and the accuracy of the results 
has been compared. The A(^f ) function can be fit in the region from to 0.5 
with a quadratic form with good accuracy, 

(4vr) 4 A(0 < x < 0.5) w -0.005648 - 0.000525x + 0.005195x 2 . (42) 
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B Appendix B 

This appendix is a minor modification of Refs. |5| relevant to our case. For 
simplicity we have taken in this appendix all the pion masses normalized to unity. 
The multiple space-time integrals will be denoted by 



(...) 



d d ri 
i(2ix) 



«-»= 'i&/s&(->- < 43 > 



We have also defined the measure as 



cwrp/2) . 

L J r(3/2 + w) v 4 ; H ' v ; 



with 



n«o= ttzL^t , P = f~l ( 45 ) 



and 

limfdd = -^da. (46) 
w _>o L J 16vr 2 v ; 

To show the procedure we will use the vertex diagram in the scalar case: Given 

the integral: 



with 





V 


= «n^» 








i=i 




p 1 = l 


-r 2 
' 1 


p 2 = i- 


(ri - q) 2 


p 3 = l 


-A 


p 4 = i- 


(ri + r 2 - 


p 2 = 1 


q 2 


= 0. 





(47) 



(48) 

We have integrated over one (r 2 ) of the two variables with one-loop techniques. 
We subtract the non-local infinities using: 

m=m+m = lj^ w -^- w - ] (49) 
The finite piece is contained in B(s) and we use a Cauchy representation for it: 

m= r-M±- s ( 5o) 

J4 a (a — S) 



This yields to: 



V= rM( —± } 1 s = ( P -r 1 ) 2 (51) 

J4 a Fi^w — s) 
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Now we introduce a Feynman parameterization in the way: 

1 /■! fXl fXn-2 



a\ . . .a n 



(n — 1)! / dx\ f dx 2 ■ ■ ■ I 
Jo Jo Jo 



1 



[aix„_i + a 2 (x n - 2 - x n -x) + • • • + a n (l - xi)] 1 
After this, one is able to write the following expression: 



(52) 



V = 2 [ [ dxi [ dx 2 X 2 (- r 1—^ — m-m ) 

Ji a Jo Jo \z—{r-\—RY\ 6 



[z-{n-Rf 

z = cr(l - x 2 ) + x\ + 2p ■ qx 2 (l - x 2 )(l - X\) 

R = -p(l - x 2 ) - qx 2 (l - Xi). (53) 

With the remaining integrals we have used the following notation: 

\z-P] m) = Fm[z] ' 



IJJJp g^g al3 + cycl. 

( [J^F ~ 4(m-l)(m-2) Fm - 2[zl - (54) 

Where the function F m [z] is defined by: 

F m [z\ = ^-T( W ) r(m : 2 .- m) ,m> 1 (55) 

1 (m) 

So the resulting integral can be expressed in terms of: 

V m [P- s] = 2 i^l f 1 dx x f 1 dx 2 x 2 P(x 1 , x 2 )F m [z]-m = 1, 2, 3. (56) 
Ji a Jo Jo 

where P(xi, x 2 ) is a xi,x 2 polynomial. The V m integrals can be split in two kinds: 

1. V 3 [P; s] is a convergent integral, and can be evaluated using some gaussian 
integration subroutine. 

2. Vi[P] s] and V 2 [P; s] are divergent. We then treat separately. 

By partial integration over x\ in ([56]) expression we deal with the recursion rela- 
tion: 

P(x 1 ,x 2 )V m [z\ = P v (l, x 2 )V m [z] - 2p ■ qmx 2 (l - x 2 )P v (x 1 ,x 2 )F m+1 [z] 

P v (x u x 2 ) = / dyP(y,x 2 ) 
Jo 

z = a(l -x 2 ) +x\. (57) 

We use it to obtain the maximum number of finite terms. 
With the rest we have two ways to proceed: 
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1. We introduce dimensional regularization and expand 'z w , with this, it is an 
easy task to pick up the finite and infinite pieces coming from V\\z\ and 
V 2 [z]. For example in the V 2 [z] case we obtain: 

P(x u x 2 )V 2 [z\ =2C{w)T(-w) f 1 dx 2 x 2 P{l,x 2 )z w (58) 

Ji a Jo 

working the piece coming from [da] and keeping the first term in the ~z w 
expansion one can obtain the result in terms of the /3(n,m) function: 

P(x u x 2 )V 2 [z] = 2C(w) T l~ w)T [ 3 l% w p(-2w,w + 3/2) t dx 2 x 2 P(l,x 2 ) 

1 (w + 6/2) Jo 

(59) 

As one can see those functions contains double poles, that together with 
the non-local divergences (with double poles) that appear in the subtrac- 
tion procedure should cancel. The finite part is evaluated once more in a 
numerical way. 

2. This other possibility is a more aimed one. 

Consider the divergent integrals containing Vi[P;s] and l^[P;s]. We have 
defined the function: 

E(m, n)= r 1^1 f 1 dx(l - x) m F n [z], n = 1, 2. (60) 
J4 a Jo 

By partial integration over the x variable we find the recursion relation: 

(3 + uo + m-n)E(m,n) = %-^- Q( u + 1 - n) (61) 

r(-w)r(n) 

-n{E(m, n + 1) - E(m + 2, n + 1)} 
where we have defined 

Q(a) = C 2 (^)r(-^)2 v ^4" r( r( 1 i/ ~ 2 ^~ ) a) (62) 

The function E(m,3) is finite at d — 4, so we use the previous relation to 
express E(m,n)(n = 1,2) through the divergent quantities Q and the con- 
vergent ones E{m, 3). 

After all one deals with a series in uj which looks likes: 

E(i, k) = C 2 (cu)T 2 (-uj){q(i, k, 0) + ujq(i, k, 1) + u 2 q(i, k, 2) + . . .} (63) 

A full table with the q values can be found in [P^j , we have checked their 
values and added some new ones needed for the evaluation of our integrals. 
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